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' It is a great pleasure to present this work in honor of Pran Nath, whose many 

-, ' important contributions to quantum field theory, supergravity grand unification, 

■ '^"'^ their consequences have been on the forefront of modern particle physics. 

T— I ! 

^ ' We examine the gravitational forces in a brane-world scenario felt by point parti- 

00 ' cles on two 3-branes bounding a 5-dimensional AdS space with / Z2 symmetry. 

C ^ ^ The particles are treated as perturbations on the vacuum metric and coordinate 

CO 1 conditions are chosen so that no brane bending effects occur. We make an ADM 

' type decomposition of the metric tensor and solve Einstein's equations to linear 

C " 3 ' order in the static limit. While no stabilization mechanism is assumed, all the 5D 

^ Einstein equations are solved and are seen to have a consistent solution. We find 

C " 3 , that Newton's law is reproduced on the Planck brane at the origin while particles 

^5 ' on the TeV brane a distance j/2 from the origin experience an attractive force that 

' has a growing exponential dependence on the brane position. 

Oh: 

i 1. Introduction 

_ ^ ' Models of the universe consisting of a five-dimensional bulk space bounded 

. by branes have been the subject of a great deal of study over the last sev- 

^ I eral years. The theoretical motivation for such work can be found in the 

■ - - ' Horava-Witten M-Theory (HW) i'^''^'*, which is the strong coupling limit 

of the EgxEg heterotic string theory. In HW, eleven dimensional super- 
gravity resides in the bulk that is bounded by two ten-dimensional planes, 
each of which is located at a fixed point of an /Z2 orbifold and carries 
one of the Eg, gauge groups. Upon compactification of six dimensions on 
a Calabi-Yau threefold there is an intermediate energy range where the 
universe will appear five dimensional. This is due to the fact that the inter- 
val separating the branes is 0(10) times larger than the compactification 
scale'^. The phenomenological studies of five-dimensional models were ini- 
tiated by Binetruy et.al.^ as well as Randall and Sundrum (RS)^. In these 
and subsequent studies it was shown that one can reproduce the FRW cos- 
mology at late times as well as using the location of the brane upon which 
our universe resides in the five-dimensional geometry as a possible solution 
to the hierarchy problem. 

In this paper we would like to determine the gravitational forces felt by 
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point particles placed on the branes. We will find that we will be able to 
recover the standard Newton's law only on the brane at the origin and that 
particles on the distant brane will feel a large attractive force which scales 
exponentially with its distance from the origin. 

We consider a 5D AdS bulk space bounded by two 4D orbifold planes. 
This could be a 5D reduction of Horava-Witten theory with /Z2 symme- 
try, or a two brane Randall-Sundrum model ^. We will study the specific 
case where the 4D planes have vanishing cosmological constant. The solu- 
tion of the vacuum equations has been given in ^ as 

rfs^ = e-''^^y\ijdx'dx^ + dy'' ; A{y) = (1) 

Point particles are placed on the branes and we ask for the grav- 
itational forces between them. There is a large literature on this 
subject^'^'^'^"'^^'^^'^^'^*, on the question of recovering Newton's 4D law and 
the corrections to it arising from the presence of an additional dimension. 

We assume that matter is added on the branes as a perturbation to the 
vacuum metric 

ds^ = e^'^^y {riij + hij) dx'dx^ + h^^dydx^ + (1 + ^55) dy'^ (2) 

and then solve the Einstein equations to linear order in h^v where ji = 
0,1,2,3,5, i = 0,1,2,3, and y = x^. 

The diffeomorphisms of a 5D theory with jZ^ symmetry are those of 
R\xS^ which commute with Z2. This means that for the transformation 

x^ ^x^^i^'= x'*" (3) 

one has that ^ vanishes at the orbifold points, y\ and 2/2- 

e{x\y{) = ^ = ir,{x\y2) (4) 

If one were to make a coordinate transformation with a non- vanishing ^ , 
then the branes become bent and this would create a complication when 
one imposes the boundary constraint on the brane leading to so-called 
brane bending effects. 

In previous analyses, the 5D Einstein equations are solved in Gaussian 
normal coordinates which are described by 

V = (5) 

&K,=Q^if^Kj (6) 

In general, these cannot be achieved without brane bending occurring. We 
give here an alternate analysis which avoids these complications. 
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2. Coordinate Conditions 

In general, any symmetric array hij = hji can be decomposed according to 
the ADM prescription'^ {hij = djhi) 

hij = hf^ + hfj + hij + hj,i (7) 

where hjj^ is transverse and traceless, hfj is transverse and can possess a 
trace, 

d^h^^ ^ ^ v^^h^^ (8) 
d% = ; v''hfj = f^Q (9) 



and 



One can write 



where 



hi = hj + ^hfi ; d'hj = (10) 



Oij - ^ (12) 
One can express each of these subpieces in terms of hij, i.e. 

F = ^''hij (13) 

hj/' = -Kik-Kjlh^^ - ^-Kij-Kklh^^ (14) 

h% = OikOjih'^' (15) 

hfj = O^hik - OikOjih^^ (16) 

In order to simplify the anaysis one can choose the coordinate condition 

/i5^ = (17) 

Under a linearized coordinate transformation 

a;''^a;''+e(x;y) (18) 



"■This decomposition was first introduced in (Ref.^® is a more accessible recent reprint 
summarizing the ADM formalism.) The generalization to 4-space is trivial except for the 
ambiguity in defining l/D^ in Minkowski space. However, we will always be considering 
the static Newtonian limit here where — + V^, though the size (and correct definition 
of) the higher order dynamical effect are also of interest. 
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one can always achieve h^i = using only the maintaining Eq.(4). In 
terms of u)'^ = e'^^^'^^^^{x;y) the remaining coordinate freedom is then 

co'[{x;y) = F^{x) (19) 
u;5ix;y) = -uj^5ix,y) (20) 
tJ5(.x; yi) = = ^5(2;; ^2) (21) 

where = uif + uifi, uij = OijUi^ , and Lof = u>i — Oijiv^ . is an 
arbitrary function of x (independent of y). This then allows a remaining 
gauge freedom on the metric of {A' = dA/dy): 

Shnn = -2(e-2V^),5 (22) 

Sf^ = 6A'e-^^io^5 (23) 

^(□^/i^) = 20^^^ + 2A'e-^^u^Q (24) 

ShJ =u;f{x) ; ShJ/^ = (25) 

Since by Eq.(4) we require L05{x;ya) = 0, f^{ya) is gauge invariant, and 
we can choose (using 075) 

d5f{x,ya,)^0 ; a = 1,2 (26) 

Note that since ^155 ^ 0, the invariant distance between the branes 
may change, so the equations should still be consistent even though a 
Goldberger-Wise stabilization mechanism has yet to be included in the 

analysis. However, a satisfactory RS model requires a stabilization mech- 
anism and we will see what problems the lack of such a mechanism may 
produce. 

3. Solution of the Einstein equations 

The 5D Einstein equations read 

RiJ = ^{Tii - \rhiT) ; T = r^^^T^^ (27) 

where M5 is the 5D Planck mass. Upon projecting out the transverse 
traceless modes, one is left with 

[-dl + 4A'a5 - e^^U-]hl^ = ^ E ^^'5(2/ - 2/"). (28) 

^ c«=l,2 

subject to the Z2 boundary conditions on the brane 

2A 

d.hf/' = {-ir'—^Tfj^{y^) ; a = 1,2. (29) 
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We solve these equations by Fourier transforming the a;' dependence 

y) = J dW^'fiJ^ip; y). (30) 
For the static forces on point particles on the brancs one needs (using Eq.7) 



On the Planck brane {yi 
2 



lTT 



0) we find 



hlo{j>]yi 
where 



3/3M| 



Too{yi) H Too(y2) 



D 



D 

Mil) Ji(6) 

.2 _ (p0)2 



-P 



Ni2 = 



_ ^2 (a) 



6-/i(a) 
_ Mil) 
6^1(6) 



^2(6) ^1(6) 



Mil) Mi2). 

N2{ii) iVi(6) 



(31) 

(32) 

(33) 
(34) 
(35) 
(36) 

(37) 



Mil) Mil). 

with Jfe and Nk being the Bessel and Neumann functions of order k. On 
the TeV brane (2/2 = t^p) we find 



hlo{p]y2) = - 
where 



3/3M| 



D 



D 



N2iiii,i2) = m2ii2,ii) A^22(a,6) = A^ii(6,ei) 



(38) 



(39) 



We use the other Einstein equations to determine the remaining com- 
ponents of the metric. Thus for R^i one has 

R5i = -^^^5f - 5A'dih55 + a^dshf = (40) 
The T part of Eq.(40) yields 

d5hf = 0. (41) 

Therefore hf is independent of y and can be set to zero by a gauge choice 
of (using Eq.(25) and the gauge freedom of Eq.(19)). The remaining 
terms in Eqn.(40) then give /155 in terms of 

1 



— 



3A 



(42) 
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The gauge choice of Eq.(26) then imphes the vanishing of /155 on the branes 
as weU. However, /i55 will not in general vanish in the bulk. 

The two remaining variables, /i^ and f'^ can be obtained from the R55 
and rf^Rij equations. These equations read 

{\dl AA'd,){U^h^ - \f) + e'^U-f - t^n^d,f (43) 



g2A 



and 



a=l,2 



{\dl - A'd,){U^h^ - + ^U^r - (44) 



g2A 



° a=l,2 

with boundary conditions 



X] Too{ya)S{y - Va) 



^ TooiVa). (45) 



One can use these to eliminate in terms of / 



T 



1 1 ry p2A(y') 





and the boundary conditions then determine 

tf/^(x\y„) = (-l)"Aroo(y„). (47) 

One should note that while is gauge variant in the bulk, it is gauge in- 
variant on the branes (which is why its values on the branes arc determined 
by the physical quantities Too(t/a))- We note also that the other Einstein 
equations contain no further constraints on the metric. 



4. Poles of hoo(p'^,ya) 
In Fourier space 

five) = ^^^^oM ; = V (48) 

which in the static limit gives a contribution to the 1/r Newtonian po- 
tential. Note /"^ at Ua sees only Too(j/a)) i-e-, it does not see matter on 
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the other brane. On the other hand, h^^ does see matter on both branes 
(Eqs. (32,38)) and therefore the total /loo of Eq.(31) does as weU. 

The poles of h^^f give contributions not only to the 1/r term but also 
give corrections to the Newtonian potential (these corrections 

are from the KK modes). One can see that the 1/r term in h^Q^ comes from 
the fact that as 0, 



N,^^^^-L (49) 



in the iVjj functions while the zeros of the denominator D give the higher 
order correction terms. For example in the regime ^i<d and ^2 ^1 the 
poles occur at 

^ = [(n+^)^ + ^(^)V^'^^^ (50) 



and for and ^2 ^1 the poles are at 

'^^nwe-'^y^ (51) 

In each regime, for the respective cases of the RS model and HW theory, 
the poles begin at 

(RS) mi « 7r/3e"^^^ « (10i^)(10-i'5^ = ITeV (52) 
(HW) TOi « lO^^e-'^'^f « lO^GeF (53) 

Here we have used the fact that in RS e~^^^ « 10~^^ to account for the 
gauge hierarchy. 

In order to calculate the potential we consider Kqq {m^^ya) with com- 
plex variable rn? and define 

z — 

Then by integrating g(z) over a large circle in the complex plane we can 
solve for /iQg^(m) as the sum (over m„) of the residues of the poles of Kqq . 
This results in 

h'Jivi) = -^^^[TM + e-^Py^Ton{y2)] (55) 
+ M| Ul(6)-J^(6)J "^^''^'^^ + ^1(6) ""''^'^ 
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and 

/^Sr(y2) = -^^^[Tooiyi) + e-'0y^Too{y2)] (56) 

+Mi \ Uf(6)-J?(a)J f"^^"^"^^ + ^Alar"^^"^^^ 

where the |TO=m„ indicates that the Bessel and Neumann functions are to 
be evaluated at m = m„. One can then convert the sum over m„ to an 
integral 

Y.-\dra^ — . (57) 
when the poles become dense as in the examples of Eqs. (52,53). 
5. Newtonian potential 

In the static approximation, ^ -^/p^- After a Fourier transformation, 

the part of hoo{p, Ua) then gives rise to the 1/r Newtonian potential 

in the static limit. Since /iqo = - f'^/S we can combine Eq.(48) and 
Eq.(55) to obtain for the 1/r part on the yi brane 

hooiyi) = -^j^lTooiyi) + e-'^y^Tooiy2)] + ^j^T^. (58) 

where the first term is h^Q and the second is -/^/3. Combining terms one 
gets 

^°'^y'^ = -j^^^'^y'^ - -rni^^'^^y^^ ^^^^ 

The stress tensor for a point particle is 

rjnfiu 1 ^^matter ^4/3j/2 ^^matter (60) 

and hence 100(2/) = e'^^yfhS^{x'^ — x'^) where to = TOoC^''-'' is the mass of 
a particle on the brane at position y. One then finds for the interaction 
energy between two particles 

/3mim'i 2/3mim2 

^^y'='^=-2i^-uii^ ^''^ 

where toi, fh'^ arc masses on the yi ~ brane and TO2 = moe~^y^ is a mass 
on the j/2 brane. We can then read ofl^ Newton's constant for each term 



- — between matter on the yi = brane 



I ( sttM^ ) between matter at j/i = and 2/2 
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We see that the Newtonian constant depends on whether both particles 
are on the yi = (Planck) brane or one is on the 1/2 (TeV) brane. The 
extra 4/3 factor for the 2/2 particle arises because the part of Eq.(58) 
has no Too (2/2) term. We can similarly calculate /ioo(y2) which then gives 
the potential at 2/2 to be 

where 7712 = rn2oe~^'^^ , Th'2 = mjQe"''*'^ arc masses on the j/2 brane, and 
mi = mio on the 2/1 = brane. (Here the subscript indicates the mass 
before it is rescaled by the exponential factor.) In V{y2) the e^^^^ in the 
second term (arising from -/^/3) is not absorbed in rescaling the masses 
on the TcV (1/2) brane in the second term which means that this term will 
dominate the potential giving an enormous attractive force that differs from 
the standard Newtonian form. 



6. Conclusions 

We have examined gravitational forces in a 5D model with two 3-branes 
and / Z2 symmetry. The coordinate conditions we use, hr-,^ = 0, produce 
no brane bending effects. While a stabilizing field has not been included in 
the analysis, the analysis is still consistent since /155 ^ and so the branes 
are not constrained to be fixed. (All the Einstein equations are seen to give 
consistent solutions.) 

The static Newtonian potential comes from two parts of the 4D piece 
of the metric: hj^ and While hj^ sees matter on both branes, at 
ya sees matter only on the y^ brane. In RS, the f'^ contribution on the 
Planck brane is what is needed to give the correct relation between the 5D 
Planck mass, M5, and the Newton constant Gn, i-e. 

for forces between Planck brane particles. However, on the TeV brane, 

even after the masses are appropriately rescaled, the f'^ contribution on the 
TeV brane gives an attractive force proportional to e^^^^ » 1. It remains 
to be seen if this problem can be eliminated by including a stabilizing 
mechanism, e.g., a Goldberger-Wise scalar field in the bulk which modifies 
the vacuum metric. This would then modify the perturbation to the metric 
/loo produced by matter on the branes. 
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